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ABSTRACT 

We present an analysis of a large sample of AGN environments at z 1 using stacked 
Spitzer data at 3.6 /xm. The sample contains type-1 and type-2 AGN in the form 
of quasars and radio galaxies, and spans a large range in both optical and radio 
luminosity. We find, on average, that 2 to 3 massive galaxies containing a substantial 
evolved stellar population lie within a 200 — 300 kpc radius of the AGN, constituting 
a > 8-cr excess relative to the field. Secondly, we find evidence for the environmental 
source density to increase with the radio luminosity of AGN, but not with black-hole 
mass. This is shown first by dividing the AGN into their classical AGN types, where 
we see more significant over-densities in the fields of the radio-loud AGN. If instead 
we dispense with the classical AGN definitions, we find that the source over-density 
as a function of radio luminosity for all our AGN exhibits a positive correlation. 
One interpretation of this result is that the Mpc-scale environment is in some way 
influencing the radio emission that we observe from AGN. This could be explained by 
the confinement of radio jets in dense environments leading to enhanced radio emission 
or, alternatively, may be linked to more rapid black- hole spin brought on by galaxy 
mergers. 
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1 INTRODUCTION 

It is now widely accepted that high-luminosity Active Galac- 
tic Nuclei (AGN) harbour accreting super-massive black 
holes implying that their host galaxies are amongst the most 
massive objects in existence at their respective epochs. In- 
deed, many studies have now shown that AGN preferen- 
tially res ide within fields con t aining over-densi t ies of galax- 



ies ( e. g. iHall fc GreenI \l99i: iBe st et al.l l2003|; IWold et all 
l2003l : iHutchings. Scholz fc Biaiic hi 2009). Together these 
points support the idea that AGN can be utilized as sign- 
posts to extreme regions of the dark matter density and thus 
the most massive dark matter haloes (e.g. [ Pentericci et al.l 
I2OOOI : llvison et al.1 12OO0I : IStevens et all |2003| ) at any given 
epoch. Combining this technique with large multiwave- 
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length surveys, like the Sloan Digital Sky Survey (SDSS; 
lAdelman-McCarthv et al.ll2006l ) which has identified up to 
100000 broad-line quasi-stellar objects (hereafte r quasars) 
up _to the highest measured redshifts (i.e. z=6.4, iFan et al.l 
120031 . has opened up a new era in AGN research. 

Many authors have addressed the question of whether 
the environments of radio- loud AGN, such as radio-loud 
quasars (RLQs) or radio galaxies (RGs), are any different 
from those of radio-quiet AGN, such as radio-quiet quasars 
(RQQs), with conflicting results. The first work that com- 
pared directly the environments of RLQs and RQQs was 
lYee fc GreenI (| 19841 ') . in which a marginally more significant 
over-density was detected around the RLQs in their sam- 
ple of objects at 0.05 < z < 0.55. However, a later im- 
proved study with more data and refined techniques r emoved 
the significance of this result (|Yee fc GreenI 19871). More 
work on the topic was conducted bv lEllingson. Yee fc Green 
l|l99ll ) who added more faint RQQs to the lYee fc GreenI 
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l|l987h sample. As a result they reported a significant dif- 
ference in the environments preferred by RLQs and RQQs, 
with RQQs in general preferring poorer environments at 
the 99 per cen t conf idence level. At 0.9 < z < 4.2 
iHutchings et all (|l999l ) found that RLQs occupied more 
dense envi ronments in t he ne ar-in frared than RQQs. In con - 
trast, both lWold et all (|200lh and lMcLure fc Dunlod (|200ll ) 
found the environments of RLQs and RQQs to be indis- 
tinguishable at z^O.2 and at 0.5 < z < 0. 8 respe ctively. 
More recently iKauffmann. Heckman fc BestI (|2008l ) found 
that radio-loud AGN reside in environments a factor of 2 — 3 
more dense than radio-quiet AGN in a large matched sam- 
ple of SDSS emission-line AGN in the local universe. These 
results present us with a very mixed picture. However, it is 
probable that many of them suffer in some way from small 
number statistics and/or significant selection effects. Fur- 
thermore, it may be crucial to understand how the AGN 
and their environments are linked at all reds hifts given t hat 
large-scal e radio-jet activity could enhance (jWiita [2004 ) or 
truncate l|Rawlings fc JarvisI liooj ) star formation on the 
Mpc scale, and hence be a crucial ingredient in semi-analytic 
models (e.g. Croton et al. 2006; Bower ct al. 2006). 

It is still not clear exactly what triggers the radio emis- 
sion we see in radio- loud AGN, and why we do not see 
similar emission i n other wise comparable radio-quiet AGN. 
iMcLure fc JarvisI l|2004 ) showed that, on average, RLQs 
have 45 per cent more massive black holes than RQQs, 
and that most RLQs have a black-hole mass greater than 
10* M0 . These results suggest that in order to be radio-loud 
an AGN requires a certain mass of black hole. However, 
at any given black-hole mass the range of radio luminosi- 
ties spans several orders of magnitude and there are RQQs 
which have black holes equally as massive as the RLQs. Thus 
black-hole mass cannot be the only factor that determines 
radio power; perhaps the environment of the AGN also con- 
tributes. 

It is also still not fully understood w hether there 
is a true AGN radio power dichotomy (e.g. Ivezic et"al] 
20021: iLacv et al.ll200ll: and more r ecently IWhite et al. I l2007l : 
Zamfir. Sulentic fc Mar ziani 2008), and whether the differ- 
ent types of AGN are intrinsically different or whether we 
just happen to be observing some whilst they are going 
through a period of radio-loudness. Studying the cluster- 
scale environments of these objects provides a good means 
of investigating this issue since, if they are identical objects 
going through diff'erent phases then they should have sim- 
ilar environments on these scales. Alternatively, it may be 
the case that the environment is in some way linked to the 
differences that we see in the radio properties of AGN. 

In this paper we concentrate on a number density anal- 
ysis of a sample of 173 AGN fields at the single cosmic 
epoch of z ~ 1, splitting the AGN into their classical types 
(RLQs, RQQs and RGs), as well as looking for trends with 
radio luminosity and black hole mass. We use 3.6 /xm ob- 
servations allowing us to sample the peak of the rest-frame 
stellar spectrum thus maximising our sensitivity to stellar 
mass. This work presents an analysis of the environments 
of the largest, most uniformly selected sample of luminous 
AGN yet assembled at high redshift. In so doing it forms 
an extension to previous studies of the environments of 
AGN with lower radio luminosities and at lower redshifts 
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Figure 1. Redshift vs optical absolute magnitude (SDSS i band) 
for quasars from the fi fth data release of the SDSS quasar survey 
jSchneider et al.ll2005^ . The quasars in bold are those used in our 
sample, clearly showing that we span the 5 magnitude range in 
optical luminosity available at 2 ~ 1. 



(z < 0.3) performed with SPSS da ta (e.g. iBest et alllioosl : 
IKauffmann. Heckman fc Bestlliooi '). 

In Section 2 we give details of our data, in Section 3 
we discuss the source extraction, in Section 4 we explain 
the method of analysis, then in Section 5 we present our 
results followed by a discussion in Section 6 and a summary 
of the main conclusions in Section 7. Throughout this paper 
we have assumed a flat cosmology with Ho = 72 km s~^ 
Mpc"^, = 0.3 and f^A = 0.7. All magnitudes are quoted 
in the AB system. 



2 DATA 

The data presented in this paper consist of infrared images of 
173 AGN taken at 3.6 nm with the IRAC camera on-board 
the Spitzer Space Telescope. The sample is split into three 
sub-samples, all at the single cosmic epoch of 0.9 < z < 1.1: 
75 RLQs, 71 RQQs and 27 RGs. This redshift was chosen 
as it is the minimum at which there is a large enough popu- 
lation of high-luminosity quasars to allow comparison with 
the bright quasars that are observed at higher redshifts. At 
this redshift the SDSS allows us to sample over 5 magni- 
tudes in quasar optical luminosity (see Fig. [T]). This sam- 
ple thus enables us to decouple luminosity generated effects 
from evolutionary ones, something which has plagued many 
other flux density limited studies in this area. In addition, 
the targets were chosen to have redshifts optimized, within 
the chosen range for follow-up CO surveys with interferome- 
ters such as the Atacama Large Millimeter Array (ALMA). 
Observing both unobscured (type-1) AGN, in the form of 
quasars, and obscured (type-2) AG N, the RGs, allo ws us to 
test AGN unification schemes fe.g. lAntonucciiri993l ). 

Full details of the quasars will be presented elsewhere 
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Table 1. The radio galaxies used in this paper. Column 2 gives 
the observed-frame 325 MHz flux density, column 3 the spectral 
index and column 4 the redshift. The 325 MHz flux densities and 
spectral indices {Su oc u~°') are calculated by fitting a power- 
law through available flux density measurements taken from the 
NASA Extragalactic Database (NED), except for 6C* and TOOT 
objects which are not listed in NED, see text for details. The 
errors associated with the flux densities are typically a few per 
cent. 



Name 


S325 (Jy) 


a 


z 


30 175.1 


6.939 


0.85 


0.920 


30 184 


9.097 


0.87 


0.994 


30 22 


8.348 


0.90 


0.936 


30 268.1 


15.615 


0.58 


0.970 


30 280 


16.025 


0.81 


0.996 


30 289 


8.278 


0.84 


0.967 


30 343 


13.413 


0.68 


0.988 


30 356 


6.820 


1.04 


1.079 


60 E0943-I-3958 


1.182 


0.85 


1.035 


60 ElOll-l-3632 


1.190 


0.79 


1.042 


60 E1017-I-3712 


1.540 


1.00 


1.053 


60 E1019-I-3924 


1.690 


0.94 


0.923 


60 E1129-I-3710 


1.543 


0.89 


1.060 


60 E1212-I-3805 


1.408 


1.06 


0.95 


60 E1217-I-3645 


1.402 


0.94 


1.088 


60 E1256-I-3648 


1.760 


0.81 


1.07 


60 E1257-I-3633 


1.036 


1.08 


1.004 


60*0128-1-394 


1.322 


0.50 


0.929 


60*0133-1-486 


0.742 


1.22 


1.029 


50 6.24 


0.839 


0.77 


1.073 


50 7.17 


0.469 


0.93 


0.936 


50 7.23 


0.546 


0.78 


1.098 


50 7.242 


0.304 


0.94 


0.992 


50 7.82 


0.371 


0.93 


0.918 


TOOT1066 


0.098 


0.87 


0.926 


TOOT1140 


0.298 


0.75 


0.911 


TOOT1267 


0.282 


0.80 


0.968 



(Jarvis et al., in prep) while a list of the RGs giving their 
main properties is given in Table [1] The RGs properties 
are taken from NASA Extragalactic Database (NED), ex- 
cept for th e 6C o b jects where the redshifts ar e taken from 
iBest et al.l l| 19961 ). iRawlings et"al] l|200ll ) and llnskip et al.] 
1 2005h. and the 6C* a. nd T OOT objects which are d escribed 
by I Jarvis etHI (|200ll ) and IVardoulaki et al.1 l|2009l ) respec- 
tively. Further details of the RGs will be presented elsewhere 
(Fernandes et al., in prep). 



2.1 Selection 

The quasars were selec ted by their opt i cal co lours in the 
SDSS Quasar Survey (|Schneider et al.l 120051 ). Using the 
SDSS to select the quasars allowed us to select a large 
enough initial sample that the RLQ and RQQ samples were 
chosen in identical ways. The initial sample that met the 
SDSS colour criteria for quasars was then cr oss referenced 
with the NRAO VLA Sky Su rvey (NVSS; ICondon et al.1 
Il998l ). the VLA FIRST survey jBecker et al.lll995l) and the 
Westerbork Northern Sky Survey (WENSS: lRengelink et al] 
[liil) to pick out the RLQs and RQQs. 

RLQs were chosen to have a low frequency WENSS 
(325 MHz) flux density of greater than 18 mjy which is the 
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Figure 2. Optical apparent magnitude (SDSS i band) vs radio 
apparent magnitude (NVSS 1400 MHz) for the quasar samples. 
RLQs are plotted as diamonds while RQQs a re shown as upper 
limits. The line shows the parameter i?; = 1 ijlvezic et al.ll20o3 ) 
which is used to determine radio-loudness, i.e. objects falling 
above the line are classified as radio-loud while objects falling 
below the line are classified as radio-quiet. The plot shows that, 
by this definition, all but 4 of our RLQs would be classified as 
radio-loud and that at least two thirds of our RQQs would be 
classified as radio-quiet. 




Figure 3. Stacked and averaged radio image of the positions of 
our RQQs using data from the FIRST survey. The contours show 
the 3-, 4- and 5-cr levels. Black contours are negative and white 
positive. The image dimensions are 1.5 X 1.5 arcmin. 



5-cr limit of the survey. At z ^-^ 1 this corresponds to a ra- 
dio luminosity almost entirely within the radio- loud domain. 
Fig. [2] compares our s ample to an altern ative definition of 
radio-loudness used by llvezic et al.1 (|2002l ). Here radio-loud 
objects are defined to have Ri > 1 where Ri = log(Fradio/-Fi) 
and Fradio and Fi are flux densities measured at 1.4 GHz and 
in the i-band respectively (K-corrections are not applied). 
With the exception of 4 objects, all of our RLQs would be 
considered radio-loud using this definition. We note that 



4 J. T. Falder et al. 2009 



O Radio-quiet qui 



1 .00 

Redshift 



Figure 4. Optical absolute magnitude (SDSS i-band) vs redshift 
for the quasars. RLQs are shown with asterisks and RQQs with 
diamonds. 
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Figure 5. Low frequency radio luminosity vs redshift for our sam- 
ple. RLQs are shown with open circles and RGs with filled circles 
(data are rest-frame 325 MHz from WENSS). For the RQQs, 5- 
cr upper limits (extrapolated to rest-frame 325 MHz) from the 
FIRST survey are shown. Where WENSS data were unavailable 
for the RLQs due to sky coverage (about 10 objects) the 325 MHz 
flux density was extrapolated from the NVSS survey at 1400 MHz 
assuming a = 0.7. The dashed line shows the average 5-(t limit of 
the WENSS survey, converted to a luminosity at 2; = 1 by assum- 
ing a = 0.7; the RLQs were selected to have radio luminosities 
falling above this line. The dotted line shows the average 5-cr 
limit of the FIRST survey, extrapolated to 325 MHz and again 
converted to a luminosity; the RQQs were selected to have a radio 
luminosity falling below this line. The assumed spectral indices 
for some conversions explains why some objects fall between the 
lines on this plot. 



these four objects have only one flux density measurement 
at radio wavelengths, and hence may in fact fall above the 
line if their spectral indices differ from the value we have 
assigned to them, i.e. 0.7 which is the mean value of the 
measured spectral indices. Using a low frequency radio flux 
to define the RLQs allows them to be compared more easily 
to the RGs without a severe orientation bias. 

The RQQs were defined as being undetected by the 
FIRST survey at the 5-cr level. FIRST was used for this defi- 
nition because it provides a more sensitive flux density limit 
than WENSS. In this case, the RQQs are not selected to 
be true radio-quiet objects as defined by the radio-loudness 
parameter for example. In order to get an estimate for the 
average radio power of the RQQs in our sample we have used 
a technique of stacking the radio images of non-detections 
to reveal the aver age value of the FIRST radio power (e.g. 
|WhiteeiaD|2003). In essence it involves stacking up the 
radio images at the known positions of the RQQs, weight- 
ing each image by its standard deviation and then comput- 
ing the average radio emission or obtaining a sensitive up- 
per limit. The stacked image is shown in Fig. |3l Using this 
technique we find an average flux density for our RQQs at 
1.4 GHz of 0.10 ±0.02 mjy (i.e. a 5-cr detection). Assuming 
a spectral index of 0.7 allows us to extrapolate to a 325 MHz 
flux density of 0.30±0.06 mJy which at 2; ~ 1 corresponds to 
a 325 MHz luminosity, Logio(L325/W Hz~^ sr"^) = 23.02. 

Lists of around 75 RLQs and 71 RQQs were chosen for 
observation to be matched in optical luminosity and span 
the full 5 optical magnitudes available; the selected sources 
are shown in bold in Fig. [1] The distribution of optical 
magnitudes within the selected redshift range is shown in 
Fig. 31 The RGs were selected from the low frequency, (178 
or 151 MHz; orientation indep endent) radio sample s 3CRR 
Laing. Riley fc LongaiJ Il983l ). 6CE jEaled 1 19851), 7CRS 



Willott et al.lll988l ) and TOOT surveys (|Hill fc Rawlingj 



l2003l) . Together, these surveys give 27 RGs in the same 
0.9 < 2; < 1.1 redshift range as our quasars. The reason 
therefore for our substantially smaller RG sample is purely 
due to the limit of the known RG population at z ~ 1. 

The RLQ and RG radio luminosity distribution within 
the selected redshift range is shown in Fig. [S] which shows 
that, on average, the RGs are more radio-luminous than the 
RLQs, albeit with a significant overlap. Using the FIRST 
radio images we can also place an upper limit on the ra- 
dio emission of each RQQ; see Fig. [5] and Fig. [S] which 
shows that at least two thirds of our RQQs and most likely 
more wou ld be classi fi ed as radio-quiet using the defini- 
tion from llvezic et al.l l|2002l ). In comparison to these lim- 
its the least radio-loud RLQ has a 325 MHz luminosity of 
Logio(L325/W Hz~^ sr~^) — 24.5, showing that there is at 
least an order of magnitude difference (several between the 
means) in the radio emission of our RQQs and RLQs. The 
reason for the gap between the RLQs and RQQs in radio 
luminosity seen in Fig. [5] is due to the difference in the sur- 
vey depths of the WENSS and FIRST surveys from which 
they were selected. The 5-cr limit used for the RLQ lower 
limit and the RQQs upper limit are not identical, as shown 
by the dashed and dotted lines in Fig. [5] This leaves a small 
region on the radio luminosity axis uncovered by our data, 
but it is important to note this is due to our selection rather 
than evidence for a real radio power dichotomy. 
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2.2 Observations and data reduction 



3.1 Source cuts 



We omitted from our Spitzer observations all targets which 
lie within the overlap regions of the SDSS, the Spitzer 
Wide- Area Infrared Extragalactic Legacy Survey (SWIRE; 
iLonsdale et~al]|2003f) an d the Extragalactic First Look Sur- 
vey fXFLS: iLacv et all l2005l l regions. This was because 
[Richards et alj (|2006l ) have already compiled the data for 
these objects, of which two are RLQs and 15 are RQQs. We 
also found that two of the 3C RGs have adequate data in 
the archive: 3C 356 (ID3329; PI Stern) and 3C 184 (ID17; 
PI Fazio). These data were downloaded from the Spitzer 
archive and added to our own observations which are de- 
scribed below. 

Our IRAC observation were carried out in all four 
bands between 2006 August and 2007 August. We performed 
5pt Gaussian dithers with the medium cycling pattern and 
12 sec frame-time to ensure good scattered light rejection 
and good photometry. The data were reduced with the stan- 
dard pipeline version S15.0.5 giving final maps with a pixel 
scale of 1.2 arcsec. In this paper we use the 3.6 pim channel 
to study the environments of our AGN since it samples light 
emitted at 1.8 at z 1 which is close to the peak of the 
stellar emission from galaxies. It is also the most sensitive 
IRAC channel. An aperture correction of 1.48 was applied 
to the 3.6 flux densities, as det ermined by the Spitzer 
First Look Survey (|Lacv et al.ll2005l ). Full details of the ob- 
servations and data reduction are given in Jarvis et al. , in 
prep. 



3 SOURCE EXTRACTION 

The images were cut down using iraf tasks to leave 3.4 x 3.4 
arcmin square images centered on the position of the AGN. 
This process removes any edge effects and underexposed 
edges caused by the dither pattern, but leaves enough area to 
sample the AGN environment out to a radius of 800 kpc at 
2 ~ 1. To create catalogues o f the images we used the SEx- 
TRACTOR software package ijBertin fc Arnout j Il996t ) . We 
used a detection threshold of 3 adjacent pixels each at 1.5 a 
above the RMS background level. The SEEING_fwhm was set 
to 1.67 arcsec which the measured resolution of our images. 
On inspection of the output catalogues it was found that the 
default value for the deblending parameter (0.005) resulted 
in non-detection of faint sources close to brigh t AGN and 
stars. We therefore adopted the value used by iLacv et al.l 
l|2005h of 0.0001 which gave a significant improvement, al- 
though in a small number of cases it led to the inclusion of 
spurious sources identified with diffraction spikes. 

To check for spurious sources we inverted the maps and 
ran the same source extraction configuration. The results 
showed that 96 per cent of the fields had no spurious 5 
sigma sources. However, the remaining 4 per cent had spu- 
rious sources associated with the diffraction spikes of bright 
stars (which give pronounced negative artifacts). We there- 
fore used this method to indentify those images affected 
by diffaction spikes and then manually checked their non- 
inverted catalogues for sources that were obviously spurious, 
these were then removed from the catalogues. 



The catalogues were filtered to reduce the noise produced 
by spurious and foreground sources. Filtering criteria were 
applied identically to the AGN fields and the blank fields 
(see below). 

In order to ensure that all images were of the same 
sensitivity a cut was made of source detections below a con- 
servative 5-cr level (see Fig. O. This level is defined as the 
5-cr level of the shallowest image in the sample thus ensur- 
ing that any sources used could be detected in any image. 
The cut is especially necessary because of the use of im- 
ages from SWIRE and the archive which are deeper than 
our own. Using this method, all sources in our fields with 
apparent magnitudes fainter than 21.1 (13.1 /iJy) were ex- 
cluded from our analysis. To put this limit into context, the 
absolute magnitude of the break in the Jf- band galaxy lu- 
minosity function at z ~ 1 is K* — —23.0 (|Cirasuolo et al.l 
2009b. An eUiptical galax y at 2 ~ 1 has 7^ - 3.6 = 0.27 
( Bruzual Sz Chariot ^ 2003) giving an equivalent break in the 
3.6 fj,m luminosity function of —23.3 or an apparent magni- 
tude of 7713 6 ~ 20.8. Thus at the 5-a depth of our survey we 
are sensitive to galaxies with 3.6 ^m luminosities of ~ L* or 
greater. We are sampling star-light emitted longwards of the 
Balmer break at z ~ 1 so any companions to the AGN are 
likely to be galaxies with a substantial old stellar population. 

In order to remove stars the SExtractor CLASS_STAR 
output parameter was used. This returns a value between 1 
and where 1 is a perfectly star-like object and is very non- 
star-like. We have cut all sources that have a clas s _Star 
par ameter greater tha n 0.8, as used by iBest et al] (|2003l ) 
an d ISmith et aLll|2000l ). We also ran the same analysis using 
0.95, another commonly adopted value, and found almost 
identical results, suggesting that most objects cut at the 0.8 
level had values greater than 0.95. 



4 ANALYSIS 

Once the catalogues were created they were searched for 
sources in annuli working out from the target AGN. The 
annuli were kept to a fixed area, rather than fixed width; 
this keeps the signal-to-noise and Poisson errors of similar 
size from bin to bin and also allows for a larger number of 
annuli. The target AGN is excluded from the galaxy counts 
as this would bias the results towards there being an over- 
density in the first bin. 

To get an estimate for the average number of counts 
in the field to compare with the AGN fields we have made 
use of the extra adjacent field that comes with Spitzer im- 
ages. This is a result of the way Spitzer works. IRAC has 
four detectors but can only point two at the target at any 
one time so the telescope has to offset to allow the other 
two detectors to image the target. While this happens the 
untargeted detectors image a region of an adjacent field to 
the same depth as the target fields. As these fields are not 
targeted at the AGN they can be used as blank or control 
fields, and were thus treated in exactly the same way as the 
AGN fields. Hence the region with the same exposure time 
as the main fields was source extracted in the same man- 
ner as the AGN fields and the average source density was 
computed. 
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We can also use the blank fields to get a measure of the 
local foreground and background in each region and subtract 
the mean source density in each blank field from its AGN 
field. In this case, we are measuring the source over-density 
rather than just the source density, which better allows us to 
stack the results. One possible pitfall of this approach is that 
the blank fields might be close enough to the AGN that any 
over-density will extend into them. The proximity of blank 
fields to the target fields is therefore a trade-off between the 
desire to subtract a local foreground whilst not wanting to 
be so close that any over-density is also subtracted. 

For comparison we have calculated the global source 
density level for the SWIRE fields using the same source 
cuts described above. Our background source density is 
7.34 ±0.35 arcmin^^ which compares to the overall SWIRE 
average from the 3 northern fields (covering ~ 25 deg^) of 
6.79±0.34 arcmin"'^. The quoted uncertainties in these back- 
grounds are calculated by adding the Poisson and cosmic 
variance errors in quadrature. Cosmic variance errors were 
calculated by placing 173 tiles on to the SWIRE fields, each 
having an equal area to one of our blank fields to give a 
mock survey of equivalent area to that of our real survey. We 
repeated this procedure 43 times and calculated the mean 
background level and its standard deviation which is the cos- 
mic variance error. At face value, our fields have a slightly 
higher average background level than the SWIRE fields but 
they are entirely consistent within the errors. 

Does this analysis imply that any over-density sur- 
rounding our AGN extends into the blank fields? This is 
unlikely for the following reasons. Firstly, there is evidence 
in the literature that over-densities for the most powerful 
radio galaxies (a t z ^ 1.6) extend out to, at most, 1.6 Mpc 
l|Best et al.l2003l ) whereas our blank fields are at their closest 
point ~2.8 Mpc (at z~l) from the AGN. Secondly, when we 
break the blank fields into strips we see no increase in source 
density towards the AGN. It is possible that the increased 
background could be due to stars that are not removed by 
the CLASS_STAR parameter cut. This may be a larger effect 
in our images since the SWIRE fields are confined to high 
Galactic latitudes which are often pre-selected to have low 
foregrounds whereas our fields are more evenly spread over 
the northern SDSS regions. 

4.1 Completeness 

In order to measure and correct for the completeness of our 
data we conducted extensive simulations to add and recover 
artifical sources. We added 1000 sources in to each of our 
AGN fields for 30 fiux bins (i.e. 30000 sources per image) 
and proportionally less for the blank fields as they are of 
a smaller area. These sources were chosen to be Gaussian 
with a FWHM of 2 arcsec, and were scaled to the required 
fiux. They were added across each image in a grid pattern 
in batches of 100 for the AGN fields and proportionally less 
for the blank fields. The sources were allowed to randomly 
move around within the grid as far as possible without two 
sources ever being placed close enough as to not be de- 
blended. Sources were considered to be recovered if they 
were found in the SExtractor catalogues within 1.5 pixels 
(1.8 arcsec) of their input position and had an extracted flux 
within a factor of 2 of the input value. In the AGN fields 
we computed the completeness for each annulus individu- 
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Figure 6. The apparent magnitudes of all sources at 3.6 
vs the error on that magnitude. The deeper SWIRE images are 
labelled, and the lines show the 0.2 magnitude error (5 cr) level and 
the corresponding magnitude for the shallowest image, i.e. 21.1 
mag. The two other smaller bands shown and labelled are from the 
two RGs which had already been imaged by other programmes; 
see the text for details. 
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Figure 7. An example of a completeness curve with empirical 
model fit. Our flux limit is 13 /^Jy which is therefore the lowest 
flux at which we use a completeness correction. 



ally which allows us to correct for bright objects masking 
regions of certain annuli. In the blank fields we used the 
average completeness for the whole blank field region as we 
are only interested in attaining one value per field. 

To eliminate the scatter in the measured completeness 
curves we then fitted them with an em pirical model of the 
form Completeness = [S"-) / {b + cS"-) ijCoppin et al.l 120061 ) 
where S is the 3.6-/^m flux density and a, h and c are con- 
stants that are fitted. An example of one of these fits with 
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Figure 8. Plot showing tlie optical absouluto magnitude (SDSS 
i band) vs the completeness at our flux limit (13 /^Jy) in the first 
annulus for all the quasars in our sample. Correlation analysis 
shows a correlation at the 99.7 per cent level. 



the overlaid data points is shown in Fig. [T] This process 
then allowed us to apply a different completeness correction 
to every annulus in each image as a function of flux which 
corrects for any lost area due to bright stars or diffraction 
spikes as well as sources lost due to being in crowded regions. 

As expected, the completeness is lower than average in 
the first annulus of the quasar fields which is due to the 
bright AGN hindering the detection of faint sources (se e 
Fig. EJ. This effect was also noted bv lYee fc GreenI (Il984 ). 
Indeed, we find a highly significant correlation between the 
optical magnitude of the AGN and the completeness of the 
first bin for the quasars at the 99.7 per cent level using cor- 
relation analysis (Spearman rank and Kendall tau). Hence 
applying the completeness corrections was found to boost 
the source density in the first annulus relative to the others 
(although the significance is not changed as the error is also 
scaled by the completeness). This effect is not as prominent 
in the RG fields as they are generally less luminous in the 
optical/infrared, due to the quasar nucleus being obscured 
so that we only see the stellar light from the host galaxy. 



5 RESULTS 

Conducting the radial search on individual fields does not 
give significant results because the IRAC data are domi- 
nated by foreground sources. There are also large Poisson 
errors resulting from counting small numbers of objects. 
However, these problems are reduced if multiple fields are 
stacked up because the larger number of counts involved re- 
duces the associated Poisson errors, as well as averaging out 
the variations in foreground sources. Hence in this paper we 
concentrate mainly on a statistical analysis of stacked data 
from multiple fields. 
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Figure 9. Histogram showing the average source over-density 
in all the AGN fields after they were corrected for completeness 
and the local background has been subtracted. The error bars for 
each field arc the Poisson error for each bin combined in quadra- 
ture with the Poisson error on the blank field level that was sub- 
tracted, both scaled by their mean completeness correction; these 
values are then added in quadrature to give the final error bar on 
the stacked histogram. The dotted line simply highlights the zero 
level. 



Table 2. The over-density shown in Fig.[9]tabulated for each bin 
and as an average for the outer bins (i.e. excluding bin one). The 
table shows the over-density (AA^), the over-density error (err 
AA') and the Sigma value which is the number of 1-sigma error 
bars (err AA) that the over-density (A7V) is above zero. 



Bin 


AA 


err AN 


Sigma 




(arcmin~^) 


(arcmin~^) 




1 


1.80 


0.21 


8.50 


2 


0.53 


0.20 


2.71 


3 


0.40 


0.19 


2.08 


4 


0.28 


0.19 


1.46 


5 


0.22 


0.19 


1.17 


6 


0.45 


0.20 


2.29 


Outer bins 


0.38 


0.09 


4.26 



5.1 Over-density for the whole AGN sample 

In Fig. |9]we show the results of a radial search conducted on 
all our AGN fields, centered on the AGN, stacked together 
and averaged. Note that the local background level calcu- 
lated from the blank field of each AGN has been subtracted 
so the plot shows source over-density rather than source den- 
sity. The error bars for each field are the combined Poisson 
errors of the count in the AGN field added in quadrature 
to the Poisson error of the blank field value that was sub- 
tracted, both scaled by their mean completeness correction; 
these are then added in quadrature to give the final error 
bar on the stacked histogram. It is obvious that there is an 
over-density in the AGN fields as a whole, since all of the 
bins have a nominal value above the background level, but 
also a big increase in this excess towards the position of the 
AGN. 

It thus appears that the bulk of the over-density is con- 
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Figure 10. Histograms showing the averaged source over-density for the AGN sub-samples. The right-hand panel shows the RQQs, the 
middle panel shows the RLQs and the left-hand panel shows the RGs. The error bars for each field are the Poission errors for each bin 
combined in quadrature with the Poisson error on the blank field level, that was subtracted, both scaled by their mean completeness 
correction; these values are then added in quadrature to give the final error bar on the stacked histogram. The dotted lines simply show 
the zero levels. 



centrated within the first bin or within 300 kpc (physical 
units) of the target AGNs. The chosen bin size came from a 
trade off between gaining a reasonable signal-to-noise in the 
first bin whilst not watering down the over-density by includ- 
ing too large an area. The radial search was thus repeated 
using annuli of varying size and it was found that the opti- 
mum size is that adopted, which has an area of 1.4 arcmin^. 
This exercise also showed that, even when smaller bins were 
used, the central over-density still extended to 200 — 300 kpc 
before dropping towards the background level. 

To calculate the significance of the over-density in each 
bin, we simply use the number of 1-cr error bars that the 
over-density is above zero. This is shown in Table [5] which 
shows the relevant values for each annular bin of the his- 
togram shown in Fig. [9] It can be seen that if we just con- 
sider the first bin this is over-dense at the 8.5-0" level. 

In order to test whether the over-density in the outer 
bins (excluding bin one) is significant we added the counts in 
the outer bins and calculated the combined Poisson error. 
This gives an over-density of 0.37 ± 0.09 arcmin"'^ and a 
significance of 4.3 a (see Table [2]). Therefore, the outer bins 
are also significantly over-dense (out to at least 800 kpc) 
but the bulk of the over-density is mainly found in the first 
bin. This pattern of a sharp peak in the central source over- 
density and then an ex tended flatter over-density was also 
reported bv lBest et al.l (.20031 for po werful radio galaxies at 
z ~ 1.6, and bv lSerber et al.1 l|2006l ') for M ^ -22, z ^ 0.4 
SDSS quasars. 



5.2 Over-density versus AGN Type 

To investigate trends with AGN classification we have di- 
vided the AGN into RQQ, RLQ and RG sub-types and per- 
formed a similar analysis to that described in Section 15.11 
The results of this analysis are shown in Fig. 1101 On first 
inspection, the most notable difference between the subsets 
is the apparently larger over-density in the RLQ and RG 
samples when compared to the RQQ sample. 

Using the same method as for the whole sample de- 
scribed in Section 15.11 we have calculated the significance 
of the over-density in each annulus for each of these sub- 



Table 3. The over-densities from Fig. llOl shown for the first three 
bins and as an average for the outer bins (excluding bin 1). The 
table shows the over-density (AAf) and the sigma value which 
is the number of 1 cr errors that the source density is above the 
blank field source density. 



AGN type 


Bin 


AAf 
(arcmin~^) 


Sigma 


Radio 


1 


2.51 


4.01 


galaxies 


2 


1.73 


2.81 




3 


-0.22 


-0.41 




Outer bins 


0.52 


2.15 


Radio-loud 


1 


2.05 


5.90 


quasars 


2 


0.33 


1.05 




3 


0.74 


2.29 




Outer bins 


0.37 


2.73 


Radio-quiet 


1 


1.25 


3.67 


quasars 


2 


0.29 


0.89 




3 


0.28 


0.89 




Outer bins 


0.34 


2.48 



samples. The main results are summarized in Table|3]which 
shows the over-densities in the first three bins and their sig- 
nificances and the same values in just the outer bins (ex- 
cluding bin one). Focusing on the first bin, which in all cases 
appears to contain the bulk of the over-density, we find over- 
densities for the RGs of 4.0 cr, for the RLQs of 5.9 a and 
for the RQQs of 3.7 a. The lower significance level for the 
RG sample probably just reflects the small sample size. It is 
worth noting that in this subset the second bin is also sig- 
niflcant at the 2.8 cr level, which if combined with the first 
bin would give a more significant over-density. 

To quantify the difference we see in the first bins of 
the RLQs and RQQs we conducted a Mann- Whitney test, 
which is a non-parametric form of a t-test, and gives the 
probability that two data sets, in our case the source over- 
densities around our AGN, have the same mean. Considering 
source densities in the first bin this test indicates that the 
RLQs inhabit, on average, more dense environments than 
the RQQs at the > 96 per cent confidence level. It therefore 
appears that all of our AGN have a significant over-density 
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Figure 11. Histograms showing the distributions of RLQ and 
RQQ black-hole mass estimates. 



out to ~300 kpc with the radio-loud objects, on average, 
having larger over-densities. Interestingly the RGs appear to 
have a larger central over-density than in the quasar fields 
although more RG data are needed to confirm this tentative 
result. 

In order to test that these results are not being caused 
by the presence of a few very over-dense fields which happen 
to fall into our radio-loud samples we removed the 5 most 
over-dense fields from our RLQs and the five least over-dense 
fields from our RQQs. The observed trend for the first bin 
still holds, although obviously there is a drop in significance 
due to removal of the most extreme fields in each case. It is 
clear from this analysis, however, that there is a fair amount 
of field-to-field variation which could be attributed to either 
actual changes in the enviromental richness between fields or 
to field-to-field changes in the foreground/background con- 
tamination. 



5.3 Black- hole mass vs environmental density 

In order to try to understand the contrast between the fields 
of the radio-loud and radio-quiet objects we can use esti- 
mates of black-hole mass to search for possible trends. For 
example, are the biggest over-densities found around the 
AGN which host the largest black holes? 

Our black-hole mass estimates are computed using the 
virial estimator and the Mgll line at 2 800 Ausing SPSS 
spectr oscopy, a technique desc ribed by McLure fc Jarvid 
l|2002[ ). and based on work by iMcLure fc Dunlod (|20oi v 
Note that we do not have black-hole mass estimates for 
the radio galaxies because the broad-line region is obscured 
in these objects, but estimates are available for all of the 
quasars. 



To test whether the difference we find between the 
fields of the RLQs and RQQs might be related to intrinsi- 
cally different black-hol e mass distribut i ons, as suggested 
in pr evious papers fe.g. iLacv et al.|[200ll : iMcLure &: Jarvid 
|2004| '1. we again conducted a Mann- Whitney test. The 
test suggested that the mean black-hole masses for the 



two samples are not significantly different; the returned 
probability that they have the same mean is 0.32. Moreover, 
a two-sample Kolmogorov-Smirnov (K-S) test returns a 
probability of 0.90 that the samples are drawn from the 
same parent distribution. We stress that this result does not 
contradict previous work but is rather a direct consequence 
of the manner in which our samples were initially matched 
in absolute optical magnitude and colours as well as the 
relatively small sample size. The mean black-hole masses 
are (Logio(MBH/ Mq)) = 8.87 ± 0.06 for the RLQs and 
(Logic (Mbh/M©)) = 8.81±0.06 for the RQQ so the means 
of the two samples are consistent and well within 1-a of 
each other. In any case the real uncertainties are likely to 
be bigger once systematics are taken into account. The 
distribution of black-hole masses is shown for both samples 
in Fig. [Til 

We can go one step further by testing whether over- 
density is related to black-hole mass regardless of AGN clas- 
sification. This is shown in the left panel of Fig. [121 To search 
for possible trends we performed a correlation analysis on 
the AGN over-density measurements in the first bin. We 
have used two types of correlation analysis: the Spearman 
rank order test and the Generalized Kendall's tau test. Both 
tests suggest a low probabilty of a correlation (see Table [4}. 

A possible source of uncertainty in the black-hole 
masses of RQQs and RLQs could be their orientations 
with respec t to the observer. As sho wn bv lJarvis fc McLurd 
(i2003 ) and jjarvis fc McLurd (|2006l 'l one would expect that 
sources with bright core radio emission would be preferen- 
tially aligned pole-on to the observer. Coupled with a disc- 
like broad-line region, such a bias would result in lower de- 
rived black-hole masses for the RLQs relative to the RQQs 
given the same optical selection. However, our initial se- 
lection, based on low-frequency radio emission using the 
WENSS at 325 MHz, means that we minimize beaming ef- 
fects as the radio emission at low frequencies is dominated 
by the extended, optically thin, lobe emission. Therefore, 
although we cannot rule out completely a link between en- 
vironmental density and black-hole mass for our sample, it 
would seem very unlikely. 

5.4 Radio luminosity vs environmental density 

Having upper limits for the RQQ radio power (see Sec- 
tion I2.1|l allows us to investigate the environmental den- 
sities of all our AGN as a function of radio luminosity (see 
Fig. I12|) . Since the RGs are not selected in the same way as 
the quasars, caution must be exercised in their inclusion in 
this analysis. It is also worth noting the size of the error bars 
on the over-density in these figures, as shown by the mean 
error bar (±3 arcmin^'^). This explains the large scatter on 
any correlation. 

The censored data (e.g. the RQQ upper limits) require 
that we use survival analysis techniques for our correla- 
tion analysis. We therefore used iraf which provides three 
types of test: the Spearman rank order test, the Generalized 
Kendall's tau test and the Cox proportional hazard model. 
All three tests can handle one type of limit in the depen- 
dent variable, in our case upper limits in the radio lumi- 
nosity. The results of the correlation analysis are shown in 
Table [4l With the exception of the Cox proportional haz- 
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Figure 12. Source density witliin ~300 kpc at z~l as a function of black-hole mass for the quasars (left-hand panel) and as a function 
of radio luminosity at 325 MHz for all AGN (right-hand panel). In the right-hand panel, the open circles show the RLQs, the filled circles 
the RGs and the upper limits the RQQs; also shown in this panel is an example of the size of the error bars on the over-densities in both 
panels. The error bars would show the Poisson error for each AGN field combined in quadrature with the Poisson error on the blank 
field level, that was subtracted, both scaled by their mean completeness correction. 



Table 4. Correlation analysis on the data from Fig. 1121 i.e. source density within ~300 kpc at z ^l vs black-hole mass for the quasars 
and then vs radio luminosity for all AGN. AGN(Lradio > 25 or 26) corresponds to AGN with Logio(I/325/W Hz~^ sr~^) > 25 or 26. 
Coefficients shown are the correlation coefficient except for the Cox hazard test which gives a value. The significance shown is the 
confidence level at which the null hypothesis (i.e. no correlation) can be rejected. Survival analysis was used on the radio luminosity data 
as these contain upper limits on the RQQs radio luminosity. 





Test 


Coefficient 


Significance 


Quasars 


Spearman rank 


0.080 


0.666 


(Mbh vs over-density) 


Kendall tau 


0.054 


0.669 


All AGN 


Spearman rank 


0.160 


0.961 


(radio luminosity vs over-density) 


Kendall tau 


0.215 


0.966 




Cox hazard 


2.109 


0.854 


AGN(L,,dio> 25) 


Spearman rank 


0.226 


0.963 


(radio luminosity vs over-density) 


Kendall tau 


0.155 


0.966 


AGN(L,,dio> 26) 


Spearman rank 


0.483 


0.992 


(radio luminosity vs over-density) 


Kendall tau 


0.121 


0.999 



aid model the tests give evidence for a correlation at the 
> 95 per cent confidence level when all AGN are included. 
We also restrict our analysis to those sources that have radio 
luminosities above the traditional separation between FRI 
and FRI I radio sources at Lokio (Z/325/WHz~^ sr~^) = 25 
(see e.g. IClewlev fc Jarvid |2004| ). At this threshold there 
is also a divergence in the space-density evolution with 
redshift, with the higher luminosity radio sources tending 
to evolve more strongly than the lower lu minosity sources 
l|Clewlev fc Jarvij|2004l : ISadler et al.ll2007D . For this sample 
we again obtain correlation at the > 95 per cent level. How- 
ever, if we employ a much more conservative low-luminosity 
cut-ofF of Logio(I/325 /W Hz~^ sr~^) = 26 which ensures that 
all of our sources lie well within the FRII regime then we 
find that the correlation is significant at the > 99 per cent 
level. 

Intriguingly, in the high radio luminosity range the 
RLQs and RGs appear to show the same trend of increasing 



source over-density with radi o luminosity, as would be ex- 
pected in the unified scheme (|Barthej^l989'). This improve- 
ment in the correlation analysis results for high radio lumi- 
nosity could of course be a real effect or instead it might be 
a result of foreground/background contamination in some of 
the low-radio luminosity AGN causing them to appear more 
over-dense th a n the y actually are. However, recent work by 
iDonoso et all (I2009D finds that the environmental densities 
of RLQs and radio-loud AGN (i.e. RGs) match only for radio 
luminosities of Logio(L325/WHz~^ sr~^) > 25.4 (after con- 
verting into appropriate units). Their interpretation is that 
the unified scheme for radio-loud quasars and radio galaxies 
might only be valid for high radio luminosit ies. Finally, we 
note with interest that iDonoso et al.l (|2009l ) also report a 
similar trend to that found here of increasing environmental 
density with radio luminosity for RLQs. 
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6 DISCUSSION 

Our results indicate that, on average, AGN at 2 ~ 1 have 2 
to 3 massive {^L*) galaxies containing a substantial evolved 
stellar population in their ~ 300 kpc-scale environments, 
representing an over-density relative to the field. 

Moreover, we find evidence that the radio emission we 
observe from AGN is in some way related to the galaxy 
density in its environment. Specifically, we find that the 
RLQs and RGs occupy more dense environments than the 
RQQs. Secondly, if we ignore AGN classification and simply 
measure the over-density as a function of radio luminosity 
(Fig. I12[) we find evidence of a positive correlation between 
the two (Table [i]). 

Initially one might expect that these results could be 
due to the radio-loud objects having intrinsically larger black 
hole masses than their radio-quiet counterparts. Indeed, re- 
cent work has found evidence for this to be the case (e.g. 
iLacv et al.ll200ll : lMcLure fc JarvidbOO^ I. If we then assume 
that the largest black holes reside in the largest dark matter 
haloes, or are found closer to the centre of their haloes, they 
would thus also have the highest density environments. How- 
ever, when we compare our results to the quasar black hole 
mass estimates we find that the black-hole masses of the 
RLQs and RQQs are statistically indistinguishable. More- 
over, the over-density around the quasars shows no signifi- 
cant trend with black-hole mass (Fig. I12|l using a correlation 
analysis (Tabled. 

At this point it is worth considering whether there 
might be a systematic bias in the way our black-hole mass 
estimates are made for the radio-loud and radio-quiet sub- 
sets. For example, if the Eddington ratios are for some rea- 
son systematically lower for the radio-loud objects then they 
must host more massive black-holes to produce a given op- 
tical luminosity. One possibility is that radio-loud objects 
are powered by a radiatively inefficient accretion process 
such as Bondi accretion of the hot phase of the IGM (e.g. 
iHardcastle. Evans fc CrostonI |2007| ) while the radio-quiet 
objects are accreting cold gas in the standard manner. How- 
ever, such a possibility is easily dismissed since radiatively 
inefficient accretion processes can only explain the multi- 
wavelength properties of low-excitation radio sources which 
are almost all Fanaroff-Riley (FR) class I objects whereas 
our entire sample is made up of QSOs, with the radio-loud 
objects all having radio luminosities typical of FR II sources 
(Fig. Furthermore, given that our samples are matched 
in absolute optical magnitude and optical colours, any dif- 
ference in accretion properties must contrive to produce a 
distribution of Mgll line widths that would lead to identical 
black-hole mass distributions; this seems unlikely. 

Therefore we are led to conclude that the environments 
of the AGN are somehow affecting the differences we observe 
in their radio properties. This is one of two possible scenar- 
ios that could explain our results, the other being that the 
AGNs radio emission is infiuencing the environmental den- 
sity. However, this is a much harder scenario to envisage as 
the ~100 kpc-scale radio jets would need to infiuence galaxy 
formation on Mpc scales. 

It is known that tidal stripping would be more prevalent 
in denser environments, as more close encounters or mergers 
with other galaxies would occur. The inter-galactic medium 
(IGM) would therefore be denser in regions of higher galaxy 



density. A higher IGM density gives more material for ra- 
dio jets to work on, increasing the radio luminosity produced 
through synchrotron lo sses. This effect, k nown as jet confine- 
ment, was discussed bv lBarthel fc Arnaud (.1996. ) who use it 
to explain the unusually steep far-infrared to radio spectral 
slope of Cygnus A as boosting of the radio luminosity caused 
by a higher environmental density. The estimated enhance- 
ment in radio luminosit y for AGN in clust e rs, co mpared to 
the field, was given bv lBarthel fc ArnaudI l|l996D as ~ 1.5 
orders of magnitude which might be sufficient to explain 
our results, although a realist ic physical model is clearly re- 
quired . In the local universe iKauffmann. Heckman fc Best! 
( 200 ^ find a similar difference in matched samples of radio- 
loud and radio-quiet emission line AGN from SDSS; they 
also offer an explanation in terms of radio jets being en- 
hanced in denser environments. Our results extend this re- 
lationship to higher luminosity objects at higher redshifts. 

The idea of jet confinement may be able to explain the 
results for the RGs and RLQs. Whether this idea can be 
extended to the RQQs, which typically have radio emission 
at least an order of magnitude lower than the RLQs, is less 
clear. Classically in unified schemes, RQQs are thought of as 
the quasars without kpc-scale radio jets, implying that they 
are physically different to RLQs. However, there is some ev- 
idence in the literature to suggest that the radi o properties 
of RQ Qs and RLQs are not so very different (see lLacv et al.l 
I2OOII for example). It has been suggested though, that the 
use of the FIRST survey, which is not sensitive to extended 
radio emissio n, may in fact have been responsible fo r these 
findings (Ivcz i'c et al.ll2002l ). However. iKukula et all (|l998l ) 
detect kpc-scale radio jets around 34 per cent of a sample of 
27 RQQs. If this is representative of the whole RQQ popu- 
lation then perhaps the classical view is no longer valid and 
there is infact more of a continuum of quasar radio proper- 
ties. This certainly makes the source of the radio emission 
from AGN easier to understand, as explaining a dichotomy 
in their properties is difficult. In our sample, there is def- 
initely a gap in the radio luminosities of our quasars (see 
Fig. O in the sense that the RQQ upper limits are typically 
an order of magnitude below the least radio-loud RLQs. This 
effect is, however, intrinsic to our sample, since we have used 
different surveys which have different depths to define the 
RLQ and RQQ sub-samples. Therefore this is not sugges- 
tive of a radio power dichotomy in the quasar population 
although at the same time it is not evidence that there is 
not one. 

An alternative explanation is that the differences we 
observe in AGN radio properties are caused by the spin 
states of their black-holes. This is certainly plausible as it 
could much more readily explain a radio power dichotomy, 
if indeed it turns out that there is one. It could also 
explain our results if for some reason the spin is affected 
by the environmental density. This idea is mentioned in 
the literature and usually takes the form of black holes 
spinning faster in dense environments due to increased 
exposure to merg ers , whi ch spin up the black hole (e.g. 
IWil son fc Colbert] Il995l. 'Moderski. Sikora and Lasotal 
1 199^. IVolonteri. Sikora fc Lasota 20oJ and 



ISikora. Stawarz fc Lasota 20071). However, the one ma- 
jor problem with black hole spin as an explanation is that 
it is as yet not observationally measurable in AGN and so 
the hypothesis is not yet testable. 
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7 CONCLUSIONS 

We have conducted an analysis of the environments of a 
large sample of AGN at z ~ 1 in order to study the rela- 
tionship between AGN activity and environmental density 
at an epoch close to the peak in AGN activity. Our main 
conclusion are: 

(i) The AGN fields show, on average, an excess above 
the field of 2 to 3 massive galaxies containing a substan- 
tial evolved stellar population. Most of this over-density is 
confined within a radius of 200 — 300 kpc of the AGN al- 
though there is good evidence for a lower level over-density 
extending out to the Mpc-scale. 

(ii) We find evidence for a trend of increasing galaxy over- 
density with increasing AGN radio luminosity. Since the 
RLQs and RQQs have indistinguishable black hole mass dis- 
tributions the observed difference in environmental density 
is not a result of observing different populations of objects. 
This leads us to conclude that the radio power of an AGN 
is in some way infiuenced by the environmental density in 
which it resides. 

(iii) Our results could be explained by the boosting of ra- 
dio jets in areas of higher IGM density which are known to 
exist in galaxy-dense regions due to mergers and tidal strip- 
ping of galaxy gas. It is unclear whether this explanation 
can be extended to the RQQs in our sample. 

(iv) It is of course entirely possible that the radio prop- 
erties of AGN are not determined by a single parameter but 
instead a combination of parameters. It is clear from pre- 
vious work that for an AGN to be radio-loud it req uires 
a certain mass of black hole (jMcLure fc Jarvid |2004 ) but 
this cannot be th e only factor involved. This work and 
that of others (e.g . iKauffmann. Heckman fc" Bestll2008l and 
IDouoso et al.ll2009l ) suggests that there is some link between 
radio-loudness and the environment as well, and there have 
been several theoretical papers proposing that black hole 
spin could also be responsible. 

Future observations, such as deep radio observations 
with the Low- Frequency Array (LOFAR), will allow us to 
investigate the environments as a function of radio luminos- 
ity well down into the radio-quiet quasar regime and up to 
higher redshifts, and place firmer constraints on our conclu- 
sions based on this sample alone. 
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